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Abstract 

A reflection-asymmetric deformed oscillator potential is analysed from the classical 
and quantum mechanical point of view. The connection between occurrence of shell 
structures and classical periodic orbits is studied using the "removal of resonances 
method" in a classical analysis. In this approximation, the effective single particle 
potential becomes separable and the frequencies of the classical trajectories are easily 
determined. It turns out that the winding numbers calculated in this way are in good 
agreement with the ones found from the corresponding quantum mechanical spectrum 
using the particle number dependence of the fluctuating part of the total energy. When 
the octupole term is switched on it is found that prolate shapes are stable against chaos 
and can exhibit shells whereas spherical and oblate cases become chaotic. An attempt 
is made to explain this difference in the quantum mechanical context by looking at 
the distribution of exceptional points which results from the matrix structure of the 
respective Hamiltonians. In a similar way we analyse the modified Nilsson model and 
discuss its consequences for nuclei and metallic clusters. 
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1 Introduction 



The development of new experimental technique in recent years has considerably increased 
the accuracy of the measurement of nuclear properties at low excitation energy. As a result, 
the superdeformed rotational bands have been discovered opening new perspectives on 
new phenomena like the one of identical bands. Low lying negative-parity states, parity 
doublets, alternating parity bands with strong El - transitions in radium and nearby nuclei 
are another interesting phenomena which could be considered as a manifestation of the 
octupole deformation. Possible occurrence of softness of superdeformed nuclei having the 
prolate shapes with respect to the octupole deformation is one of the current topics in 
nuclear structure physics (see for review f^). 

Ten years ago the pioneering experiments carried out by a group of W. Knight [^] (see 
also 1^, ^) opened a new fascinating field in study of many-body system, the atomic cluster 
phenomena, an intermediate form of matter between molecules and bulk systems. Exper- 
imental results on metallic clusters concerning abundance spectra, ionization potentials, 
photoexcitation etc was immediately interpreted in terms of electronic quantized motion in 
a spherical effective potential. 

Phenomena observed in many-body systems like nuclei and clusters could be explained 
within the mean field approach based on the symmetry breaking mechanism related to 
quantized single-particle motion. The quantisation of a system of Fermions moving in a 
common potential leads to a bunching of levels in the single-particle spectrum, known as 
shells. The high level density around the Fermi level (large level bunching) corresponds to 
less stable system. When a spherical shell is only partially filled, a breaking of spherical 
symmetry, resulting in an energy gain, gives rise to a deformed equilibrium shape. The 
basic concept of the deformation in the mean field approach is the Jahn- Teller effect 0, 
the mechanism first time proposed for molecules, which leads to a spontaneous symmetry 
violation (see discussion in [§]). 

On the other hand, the existence of magic numbers corresponding to spherical and de- 
formed nuclei/clusters could be explained in terms of classical trajectories based on the 
periodic orbit theory of Gutzwiller [|10|. According to the semiclassical theory the fre- 
quencies in the level density oscillations of single-particle spectra of nuclei are determined 
by the corresponding periods of classical closed orbits. The short periodic orbits give the 
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major contribution to the gross shell structure [TT|, |T^. Depending on the particular mean 



field potential a deviation from spherical symmetry can lead to chaotic motion in the corre- 
sponding classical problem, and the shell structure of the corresponding quantum spectrum 



is affected or even destroyed depending on the degree of chaos [0, [1^, [I5[ . 

Simple interpretation of the strong peaks, which are observed in the abundance spectra 
of metallic clusters, in terms of the spectrum of one-particle orbits in a spherical potential, 
which are associated with periodic orbits of the corresponding classical problem |16| , p!7| , , 
provide us with the important key elements. The fact that the triangular orbits play an 
important role for a family of quantum states with A/ = 3 (here / is an orbital angular 
momentum, see [^) implies an importance of octupole degree of freedom for configurations 
with partially filled shells . 



In the present paper we consider a refiection-asymmetric axially deformed oscillator 
potential and focus our analysis on the distinction between orderly and chaotic motion in 
many-body systems, like nuclei and metallic clusters. One major result of the classical 
analysis demonstrates that the prolate case including octupole deformation is still quasi- 



integrable. In this way the quantum mechanical shell structure found earlier |15] is given 



a proper theoretical foundation Similarly, we analyse the modified Nilsson model 



including the /^-term and neglecting spin-orbit coupling, and the consequences for nuclei 
and clusters. Quantum mechanical calculations confirm that the /^-term which gives rise to 



chaotic behaviour in the corresponding classical problem ||2T| may interfere with the search 
of shell structure. 

2 The Model 

We investigate the classical and quantum mechanical motion in an axially symmetric de- 
formed potential which can be expanded as: 

2 

y{r, 0) = -i^r^ (1 + a2P2{cose) + agPslcos^) + ...). (1) 

In this expansion the terms proportional to the components of the deformation tensor (oj) 
are the Legendre polynomials of z-th order, respectively. In this paper we only consider terms 



up to, and including, the third order. Rewriting the potential in cylindrical coordinates [15 



3 



we arrive at the form: 

v(,,.)=—[,^+^+x^^y (2) 

where u = ojq^Ji — 02/2, h = — «2/2)/(l + 02) and A = 0:3/(2 — 0:2) We will call b the 
quadrupole deformation and A the octupole strength. For 6 > 1 (6 < 1) the potential is 
a deformed harmonic oscillator of prolate (oblate) shape coupled to an additional octupole 
deformation of strength A. 

Classically, the trajectories can be calculated using the Hamilton formalism, by 
integrating the following two equations of motion: 

dV{g,z) 



dg 

z = (3) 

oz 

Two aspects are important: First, because of the axial symmetry, the z component of the 
angular momentum (Z^) is conserved. As we will see in the next section, a finite value of 
Iz will not modify the winding numbers of the orbits. For this reason we will only consider 
Iz = 0. Second, we notice that the potential scales as V{'-yg,'-yz) = 'y'^V{g,z). This means 
that one value of the energy will produce the same periods and trajectories as any other 
energy value. 

The use of the term r^P3(cos0) instead of r^P-^{cos6) in the potential ensures a proper 
bound state problem for any value of the quadrupole deformation, provided |A| is smaller 
than some Acrit- Here, Acrit is defined to be the least value of the octupole strength for which 
the potential tends to -00 along one direction in the g — z plane. Analytic expressions for 
Acrit and the directions along which the potential opens its valleys to —00 are given elsewhere 
pO| . There is a particular value of b, namely 6 ~ 0.58 for which valleys along two directions 
[g = 0) and {g ^ OAz) will open simultaneously when A approaches its corresponding 
critical value. We expect that this case will exhibit maximal chaotic behaviour. 

3 Classical Approach 



The coupling to the octupole potential renders the problem nonintegrable |]T5[. Despite this 



it appears that in prolate situations {b > 1), the quantum mechanical spectrum exhibits 



shells |jT5|for some particular values of the octupole strength. A periodic structure which 
appears in the energy spectrum should be a reflection of one or more shortest classical orbits 
which dominate in the phase space 0. Here we present a method which reveals that the 
classical counterpart is quasi-integrable for prolate cases. It will also establish that the 
phase space structure is very close to that of a potential without octupole term but with a 
larger value of the quadrupole deformation b. 

Our classical approach is based on the "removal of resonances" method developed in 
the secular perturbation theory |22[. In the technique, the Hamiltonian written in action- 
angle coordinates, is averaged over the faster phase. Usually, prior to such an operation, 
a canonical transformation is necessary in order to remove the initial resonance from the 
unperturbed Hamiltonian (in our case the Hamiltonian without the octupole coupling). In 
the new rotating frame, one of the phases will only measure the slow variation of the variables 
about the original resonance which now becomes a fixed elliptic point. The problem is then 
treated by averaging over the remaining faster phase. In the case of a super or hyper - 
deformed potential though, there already appears to be a clear distinction between a slow 
and fast phase, therefore the canonical transformation is unnecessary. 

The complete Hamilton function written in terms of the action-angle variables of the 
unperturbed problem reads: 

J.+ ^^. + A^^'";'P"f' I (4) 



The frequencies of the two motions in z and g can be expressed as: 
dH 



ooz{J,e) 



dJz 



to f 5 sin e,(46V2 sin^ 9^ + GbJ.J, sin^ 6^ sin' 6, - 3 sin^ 6^) \ 
T 1 + ^ . . .9. . . . 9. .3 (5) 



OH 



^ / _ sin Oz sin eg^/bJz{%bJz sin 9^ + 3 Jg sin gg) \ 
^ V 2VbTz{bJz sin^ 9, + sin^ 9^)1 ' 



Obviously, the winding number of a trajectory defined as the ratio ujg/ujz is essentially 
equal to b when A is small. This means that for b sufficiently far from unity, there will always 
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be a fast and a slow moving phase, i.e. the averaging is performed over the corresponding 
fast angle. 



3.1 Prolate Potential (6 > 1) 

For prolate deformation the fast phase is dg. Averaging the Hamilton function over Og, 
we are led to a new averaged Hamiltonian H^v which is independent of Og. This means 
that the corresponding action becomes a constant of motion in this approximation: 
Jg — Je(0). Introducing the notation = 2Jg{Qi)/{muj) — ^^(0)+p^(0)/(ma;)^, the averaged 
Hamiltonian reads in the original position momentum coordinates: 



+ IT + A^^4^ T-M--,) - 37r2Fi(-, -, 2, -> 



(7) 



62 ' 27r \e 2' 2' ' 

Here, K and 2F1 are the complete elliptic function of the first kind and the hypergeommetric 
function, respectively. In this way, the problem is effectively reduced to two uncoupled one 
dimensional cases. A remainder of the actual coupling is the fact that, through ^, the z 
motion still depends on the initial conditions of the g motion. The frequency cu^ is given by 
u;z — 27r/T where 



T = / , (8) 

with Ez — E — Eg — E — ma;^^^/2 and U{z) is the effective potential of the z motion and 
is given by: 



U{z)^ 



(9) 



The winding numbers cug/cuz can be evaluated numerically, with cug — cu. It turns out that 
for a given value of the octupole strength, the frequencies are virtually independent of ^, 
and most of the trajectories in the phase space will have basically the same winding number. 
This number can be evaluated analytically at ^ = 0: 



' + ' 1. (10) 



cu 



2 V^l + A/Aerit ^/l - A/ A 

crit , 



We display in Fig.l and Fig. 2 the effective potential and the variation of uOg/uOz for a few 
values of ^, respectively. As a function of A the winding number increases monotonously. It 
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means that for a prolate case, a superposition of quadrupole and octupole deformation is 
practically equivalent to a pure but larger quadrupole deformation. Of particular interest 
is the existence of the short periodic orbits, as they become important for the quantum- 
classical correspondence. For instance, in the superdeformed case {b = 2), according to 
formula (10), the winding number becomes equal to 5 : 2 at A = 0.66Acrit and should have 
that value for most of the orbits. As a result, the quantum mechanical spectrum is expected 



to exhibit a periodic structure ||2^ of the same kind as of a pure quadrupole deformation 
with winding number 6 = 5:2. Likewise, when A = O.SAcrit, the ratio of the two frequencies 
becomes 3 : 1 and, as a consequence, a shell structure very similar to the case 6 = 3 and 
A = 0, is likely to occur. In Fig. 3 we display the actual phase space configuration {z — p^) 
for A = 0.66Acrit and A = O.SAcrit- Of importance is the large single separatrix structure 
which covers most of the surface of section. The corresponding periodic orbits lying in the 
centre of each separatrix are displayed in Fig.4. Their winding numbers are 5 : 2 and 3 : 1, 
respectively. After having tested numerically the winding number in different regions of 
the phase space, it turns out that with the exception of the innermost zones, the winding 
numbers are very close to those given by formula (9) for about 85% of the phase space. 

The expression (10) for the ratio of the winding numbers is of course valid for all values 
b > 1 for which the approximation is meaningful. It turns out that this is the case for 
b > 1.5. The significance of this statement is far reaching: it means that even for irrational 
values of b, where the unperturbed problem does not give rise to closed orbits, the switching 
on of an appropriate octupole strength will produce a situation of an unperturbed pure but 
larger quadrupole deformation. For instance, b = ^/S and A = 0.56 should essentially yield 
the situation of pure superdeformation (6 = 2). The consequence for the quantum spectrum 
will be discussed at the end of the following section. 

As mentioned in the previous section, we restrict ourselves only to zero value for the z 
component of the angular momentum. A finite value of Iz does not change the result of the 
averaging procedure. This implies that if at some value of A the quantum spectrum exhibits 
a shell like structure for one value of the angular momentum, it will exhibit the same type 
of structure for any angular momentum and therefore produce genuine shells. We return to 
this point in the next section. 
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3.2 Spherical and Oblate Potential (6 < 1) 

In these cases, the approximation procedure described above does not achieve any decouphng 
between q and z motion. For an oblate potential, the faster phase is 6z- Since the octupole 
part of the potential is an odd function of 6'^, the result of the averaging now vanishes. If 
higher order terms in the Fourier expansion of the octupole are considered, the coupling 
between the two degrees of freedom is still maintained; the attempt to obtain an analytic 
expression for the effective potential seems to stop here. The fact that the problem remains 
coupled even perturbatively indicates that chaotic behaviour may occur with increasing A. 
In the spherical case, one could in principle perform a canonical transformation enforcing one 
of the two frequencies to become faster and then averages over the faster phase. However, 
the slower frequency turns out to be a complicated function of time, in other words, the 
ratio ujg/ujz fails to remain constant even for an individual orbit. Like in the oblate case, 
the spherical potential is expected to yield chaos for sufficiently large values of A. For 
comparison, we display surfaces of sections in Fig. (5) for the spherical and oblate case 
where the onset of chaos is obvious. We mention here that the special case 6 0.58 requires 
the smallest octupole strength to become chaotic. The onset of chaos is associated with 
the disappearance of periodic (shell) structure in the corresponding quantum mechanical 
spectrum. This feature will be dealt with in the next section. 

4 Quantum Mechanical Approach 

We write the Hamilton operator in the form Hq + \Hi in a representation where Hq is 
diagonal. The basis chosen is referred to as the basis using the asymptotic quantum numbers 
n_L, Uz and A where n± = + and A = n+ — n_ Here the quantum numbers n_|_ 
and n_ are the eigenvalues of {A^y and {AJ)'^A_, respectively, where, in terms of the 
usual boson operators and a^, we use A± = {a^ =F iay)/\/2. For a fixed value of A 
this leaves two quantum numbers (refiecting the two degrees of freedom) to enumerate 
the rows and columns of the matrix problem. For A = the diagonal entries of Hq are 
thus B^^^^ = huj{n± + 1 + {uz + l/2)/b). The matrix elements of Hi are obtained from 
those of 2 ~ (4 + a,) and ~ {A+iA+Y + A^A^Y + A+A^ + (A„)t(A+)T). To get 
the matrix elements of 1/y/g'^ + z''^ in a numerically consistent way, we first calculate the 
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matrix elements Sm^n of g^ + z"^ from which the inverse square root is obtained using S^^^"^ = 
U ■ D^^/^ • U"^ where D = W ■ S ■ U is the diagonal form of the positive definite matrix 5* 
and U is the orthogonal matrix which diagonalises S. To ensure also numerically that S 
has only positive eigenvalues it is important that the matrix for z'^ is obtained by squaring 
z and not by evaluating analytically the matrix elements from (aj + a^)^; inconsistencies 
are otherwise introduced due to truncation. In this way, we also ensure that the truncated 
matrices S~^^'^ and the representation of 2z^ — 3zg'^ commute. 

4.1 Energy Levels 

In Fig. (6a) we illustrate the spectrum so obtained as a function of A for 6 = 2. The shell 
structure at about A = 0.63Acrit and A = 0.76Acrit is clearly discernible. The fact that 
the structures do not occur exactly at those values of A as calculated in section 3.1 is 
a consequence of the approximation introduced by the "removal of resonances" method. 
However, the agreement is very good, with an error of only 5%. While Fig. (6a) presents the 
spectrum for A = 0, we illustrate in Fig. (6b) the whole spectrum which is a superposition of 
all possible A-values. The shell structure is then more pronounced which is expected since, 
according to the discussion in section 3.1, the orbits will have the same winding numbers 
independent of the angular momentum; as a consequence, the quantum spectra will have 
shell structures similar to an unperturbed oscillator. To contrast with the case 6 = 2, we 
also display spectra for 6 = 1 (spherical) and 6=1/2 (oblate superdeformed) . As predicted 
by the classical analysis, there are no obvious regions where periodic structures such as 
shells would occur. The level statistics shows that such cases are actually chaotic, in that 
their nearest neighbour distribution approaches the Wigner surmise |^ . 



4.2 Total Energy Method 

To analyse the quantum spectrum we proceed in the orthodox fashion in that the total 
energy Etot{N,\) = Y^i^ is approximated by a smooth average function S{z,X) and 
the fluctuating difference 

6E{N, A) = Etot{N, A) - S{N, A) (11) 

is then subjected to further investigation. The finding of a suitable form for the average func- 
tion S{z, A) is facilitated in our case as it is well known that the leading term of Etot{N, A), 
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as a function of N, is proportional to N"^^^. We determine the five constants ao(A), . . . , a4(A) 
in S{z, A) = J2k=o o,k{X)z'^/^ by a least square fit which turns out to be perfectly satisfactory 
for all values of < A < Acrit- 

In Fig. (7a) the fluctuating part 6E{N,X) is presented as a contour plot and refers to 
6 = 2. It displays the ranges 100 < < 700 and < A < Acrit- For A = (the bottom 
horizontal line) we clearly discern the shell structure of the plain deformed oscillator. Note 
that the sharp minima (dark shadowing) occur at the positions N where a shell is closed; 
the distances are proportional to N^. When A is switched on the shell structure persists to 
a great extent; only when A approaches its critical value (top horizontal line) the structure 
begins to be washed out. There are local minima discernible at X/Xait ~ 0.76. This is a 
reflection of the enhanced shell structure discussed in the previous section. For comparison 
we also display contour plots for b = 1.5, 1, 0.58 and 0.5. As expected, the ordered formation 
of local minima along a line A = const becomes lesser pronounced and eventually disappears 
as the spherical and oblate region is approached. The greatest extent of disorder is obtained 
for b = 0.58 where the classical case is maximally chaotic. 

4.3 Magic Numbers 

The detailed shell structure is obtained from the second derivative of 5E{N, A) to which 
we turn next. For 6 = 2, we have plotted in Figs. (8) the function g{E) — SE{N + 1, A) + 
5E{N — 1, A) — 25E{N, A) versus N^^^ for a few characteristic values of A. The peaks of 
the plots represent the magic numbers which characterise the shells, and the heights of 
the peaks reflect the energy distance from one shell to the next. In Fig. (8a) the essentially 
unperturbed result (A/Acrit = 0.15) is presented for demonstration. The second row displays 
the results for the particular values of A for which the winding numbers are 5:2 (A/Acrit = 
0.63) and 3:1 (A/Acrit = 0.76), respectively. The magic numbers and the heights of the 
peaks agree well with those which are obtained from the unperturbed (A = 0) quadrupole 
deformed oscillators (third row in Fig. (8)) with b — 5/2 and 6 = 3, respectively, at least 
for N < 700, N^/^ < 8.88. The agreement extends in particular to the respective occupation 
numbers, i.e. the degeneracies; of course, the heights of the peaks do not show the same 
regularity as the corresponding unperturbed problem; nevertheless, even for the heights an 
overall agreement prevails when comparing with Figs.(8e) and (8f) where the respective 
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unperturbed structures are displayed. For higher values of we do get deviations which 
reflect upon the fact that the system is nonintegrable and cannot give complete order in all 
its results. While the agreement for lower values of was to be expected from the discussion 
in the previous section, the extent of the agreement is rather remarkable, especially for A ~ 

0. 76Acrit where an astoundingly clean shell structure reoccurs after it partially disappeared 
for a somewhat smaller value of A. It is this recurrence of shell structure which gives rise to 
the local minima in Fig. (7) as pointed out above. 

One of the questions addressed in this paper is also the possible occurrence of supershell 
effects when a quadrupole deformed harmonic oscillator is perturbed by an anharmonic 
term, i.e. the octupole deformation. Of particular interest in this respect is Fig. (8b) which 
refers to the intermediate value A = O.TAcrit where the genuinely different orbits with winding 
numbers 5:2 and 3:1 coexist. The long wave length fluctuation could well be interpreted 
as a supershell structure, even though that it is not as clearly pronounced as in a more 
transparent integrable case 0]. Yet the difference to Fig.(8d) which refers to a larger value 
of A is striking. 

4.4 Quantum- Classical Correspondence 

In Figs. (9) we display the square of the modulus of the Fourier transform of the level density, 

1. e. the function 

Fit) = \Y,e''-'\\ (12) 



n 



The spectrum is taken for 6 = 2 at A = 0.63Acrit (Fig. (9a)) and at A = 0.76Acrit (Fig(9b)), 
both spectra refer to A = only. The pronounced peaks can be directly associated with the 
periods of the classical 5:2 and 3:1 orbit, respectively, the periods obtained from Figs. (9) 
are in perfect agreement with those of the corresponding classical orbits which are found 
numerically by integrating Eqs.(3). This is a beautiful demonstration of Gutzwiller's trace 



formula |I0[ . As expected the frequencies deviate considerably from the unperturbed values, 
i.e. from the frequency associated with 6 = 2, A = 0, but also from the frequencies associated 
with h = 5/2 or 6 = 3. In units of the unperturbed value [b = 2) we flnd T5.2 = 1.2 and 
T3.1 = 1.4; the values are larger than unity in accordance with Fig. (2). We mention that 
a similar situation is encountered also in the hyperdeformed prolate case. There the two 
main orbits have the winding numbers 6 = 7:2 (occurring at A = 0.6Acrit) and 6 = 4:1 
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(occurring at A = 0.72Acrit), respectively. In a similar vein, as was emphasised towards the 
end of Section 3.1, a shell structure just like that of a pure superdeformed prolate case must 
occur for, say, b = a/3 and A < 0.56 p3|. These findings make it quite clear that, as far 



as the model is concerned, the spectrum alone cannot distinguish between an additional 
octupole deformation and an unperturbed but larger quadrupole deformation; additional 
experimental information is needed to settle this point. 

Again we stress that the high degree of order which prevails in the superdeformed and 
hyperdeformed prolate case when the octupole term is turned on, does not exist in the cor- 
responding oblate, in fact, not even in the spherical case. There, chaotic behaviour becomes 
manifest for much lesser octupole strength, which results in a complete disappearance of 
shell structure in the quantum spectrum. 

5 Exceptional Points 

The spectra displayed in Figs. (6) exhibit avoided level crossings which are related to the 
singularities |^ of the energy in the complex A plane. These singularities are known as 



exceptional points They represent the values of A where two energy levels coalesce 

when continued into the complex plane. The connection between the occurrence of avoided 
level crossings and exceptional points is similar in nature to the connection between the 
poles of a scattering function and the resonance structure of a cross section. In the same 
way as the poles of the scattering function give rise to the shape of the cross section, the 
exceptional points bring about the shape of the spectrum. Their interplay alone provides 
the mechanism for the signature of chaos in a quantum system. With the aid of exceptional 
points, criteria for quantum chaos can be found even when the classical counterpart does 
not exist [EH . 



The levels are obtained by solving the secular equation: 

det{E ~ Ho - XHi) = (13) 

To enforce coalescence of the roots of Eq. (12) the additional algebraic equation 

^det{E -Ho- AiJi) = (14) 

must be solved simultaneously. Since the Hamilton operator is irreducible with respect to 
symmetries, the fulfilment of the two equations simultaneously is generically excluded for 
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real finite A, as this would mean a genuine crossing for two levels. It has previously been 
established that a high density of exceptional points is a prerequisite for the occurrence of 
chaos in the energy spectrum Also, it is known that the statistical distribution of the 
real parts of exceptional points is close to the distribution of avoided level crossings [Q. It 
means that for the particular values of A where the density of avoided level crossings is high 
the energy spectrum generically obeys the same level statistics as that ascribed to quantum 
chaos. We have analysed the eigenvalue problem of the operator Ho + XHi beyond the Acrit 
value. Although the region A > Acrit is of no physical significance, it gives information about 
the distribution of avoided level crossings. To exclude genuine crossings between levels we 
display the A = subspace of the energy spectrum in Fig. (10). The spectrum is extended 
into the physically forbidden region for the prolate superdeformed {b = 2) potential. It 
appears that the maximum of the distribution of level repulsions occurs in the unphysical 
region for the prolate case in contrast to the spherical and oblate case (see Figs. (6c, d)). 
This can be understood since, unlike in the spherical and the oblate situation, the first 
order perturbation term of the prolate problem vanishes. As a result, the spectrum has a 
zero derivative at A = 0, and the exceptional points cannot occur near to the A = axis, in 
fact most occur beyond Acrit- In turn, the occurrence of avoided level crossings in the physical 
region for oblate and spherical potentials indicates the onset of chaos for A < Acrit- While 
this was expected from the classical behaviour, we see here the corresponding mechanism of 
the quantum mechanical matrix problem at work. It is interesting to note the connection 
between the symmetry breaking mechanism and the distribution of level repulsions: the 
system aspires to break the spherical symmetry where the density of avoided level crossings 
is high, so as to keep order and stability. The conclusion that the prolate case exhibits 
less chaos for small values of A than the spherical and oblate case, can in fact be drawn 
from the matrix structure alone, which in turn is directly related to the gross distribution 



of exceptional points p8 



6 Modified Nilsson Model 

Used with a great deal of success in the description of deformed nuclei, the modified Nilsson 
model has recently been extended to metallic clusters p5[ (see for review 0, |^). In practi- 
cal applications the Woods-Saxon potential, also used successfully for the investigation of 
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metallic clusters, has the disadvantage that it can not be solved analytically, in contrast to 
the modified oscillator potential |]26| , |27[| . The level ordering in the Woods-Saxon case falls 
between the soft-surface harmonic oscillator (HO) and the hard-surface square well. The 
same level ordering is obtained in HO by the addition of a term and reads in cylindrical 
coordinates 

2 

-.2 

Here, / is the square of the dimensionless angular momentum operator and vu is a constant. 

-,2 

To avoid a general compression of the shells produced by I alone, the average value of 

-,2 

< I >Ar is usually subtracted. However, while in nuclear structure the interest is focused 
on the lower lying levels in the single particle potential, the relevant energy range extends 
higher up for clusters. This necessitates a re-discussion of the possible usefulness in the 
application of this model to metallic clusters. For the present discussion we leave out the 

-,2 



term < I > p9| . The classical Hamilton function can then be defined: 

= ^ + ^ + S ) - ^ ((^^^ - + + f ^''^ 

Because of the axial symmetry, is conserved. Again, the zero value of Iz is sufficient for 
all relevant aspects of our analysis. The classical problem is nonintegrable and gives rise to 
chaotic motion if f/ > and 6 7^ 1. We also note that the Hamilton function scales 
as: H{ap,af) = a'^H{p,f), provided the product EU remains constant {E is the energy). 
Since the connection between U and va is given by: 

U='^vu, (17) 

the classical behaviour for EU = const will be refiected in the quantum counterpart along 
the invariant lines Evu = const. To get an idea for which energy and coupling strength 
signatures of chaos can be expected in the quantum spectrum, a brief discussion of the 
classical phase space is appropriate ||21[. In Fig. 11 a section of the z — pz plane with the 
g = plane is displayed. The phase space is noncompact as can be seen from the figure 
where the allowed regions are shaded. The parameter EU determines the structure of the 



phase space; if it is larger than the two lines z = ±1 / \JlniU intersect the ellipse. In 
this case, chaos is clearly discernible inside the ellipse for the classical motion even for small 
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deviations from unity of b. If EU < chaos occurs inside the elhpse only for b > 1.5. For 
E < the elhpse is absent and the motion appears regular. Using the expression (17) this 
translates for the quantum levels E^ into E]\fU ~ Nuj'^vu which means that, for vu = 0.07, 
chaotic behaviour is expected for > 3 if 6 = 2, while for, say, b = 1.05 its onset should 
occur only as high as > 20. Actual calculations fully confirm these expectations. The 
fact that an arbitrary large negative energy is allowed classically, is reflected also in the 
quantum spectrum; however, if the matrix problem is truncated, there is a limitation on the 
eigenvalues of / ; in any case, this part of the spectrum is virtually uncoupled to positive 
energies and has no part in avoided level crossings. 



7 Summary 

In the present paper we thoroughly investigated a simplified model to analyse the nature 
of shell effects at quadrupole/octupole axial deformations. We left out terms in the Nilsson 
model, which, albeit physically important, are prone to blur the analysis when the interest 
is focussed on the essentials, in particular a distinction between orderly and chaotic motion 
in many body systems, like nuclei and metallic clusters. Contrary to the case of spheri- 
cal potentials our Hamiltonian is nonintegrable. Using the method based on the 'removal 
of resonances' approach ||22[, we established a connection between shell structures in the 



quantum mechanical spectrum and periodic orbits in the corresponding classical problem. 
In the prolate case the classical problem can effectively be reduced to an integrable set of 
equations, since the motion in the two coordinates becomes uncoupled. This is in contrast 
to the spherical and oblate cases, where the motion becomes chaotic when the octupole term 
is switched on. In the prolate case, at particular values of the strength parameter A the 
octupole term produces a motion which resembles to a great extent that of a pure but more 
enhanced prolate oscillator. A superposition of quadrupole and octupole deformations is 
practically equivalent to a pure but larger quadrupole deformation. This provides the clas- 
sical explanation for the existence of quantum shell structure for prolate/octupole deformed 



system within the model |15, E0[. The fluctuating part of the energy has been extracted from 



the quantum mechanical single particle spectrum using the procedure which is a variant of 



the method discussed in We have found remarkable agreement between manifestations 



of shell structure for the same values of A/Acrit and the ones which give rise to large stability 
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islands in the Poincare surfaces of sections relating to the classical short orbits with winding 
numbers 2:1 (A ~ 0), 5:2 (A/Acrit = 0.63), 3:1 (A/Acrit = 0.76). In the intermediate case 
(A/Acrit = 0.7) when the orbits with winding numbers 5:2 and 3:1 coexist the long wave 
length fluctuation could be interpreted as a supershell structure. But due to the narrowness 
of these orbits the supershell structure is not as clearly pronounced as in a more transparent 
integrable case |T^. We placed our emphasis on the superdeformed case which is the most 
interesting situation when the octupole term is switched on. Shell structure is destroyed for 
smaller values of the octupole strength for lesser quadrupole deformation or even more so 
for the spherical case. Our results imply that the shell structure is supported by the su- 
perdeformed prolate/octupole in contrast to the oblate/octupole case. This conclusion was 
confirmed by the analysis of exceptional points connected with the occurrence of avoided 
level crossings in quantum spectra. It appeared that in the spherical and oblate case, the 
distribution of level repulsions, which provide the mechanism for the signature of chaos in 
a quantum system, occurs to a much greater extent in the physically relevant region than 
in the prolate case. The analysis of the modified Nilsson model without spin-orbit cou- 

-,2 

pling leads to the conclusion that using the / -term in models for metallic clusters with 
the coupling strength similar to that in nuclear physics would conflict, at least for larger 
deformations 6, with the experimental finding of shell structure. Only if the deformation h is 
very close to unity can the regular motion associated with shells prevail, provided the shell 
number is not too large. In turn, only a much smaller coupling can yield shell structure 

-,2 

for larger deformations h. Since the / -term has only a phenomenological meaning, this 
finding simply calls for caution when transferring models from nuclear to cluster physics. 
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Figure captions 



Fig.l The potential U{z) in units of ma;^/2 for a few values of A and an intermediate 
value of ^. Variation of ^ produces similar shapes. 

Fig. 2 The ratio oOp/uz as a function of A for various values of ^. A genuine ^ dependence 
is discernible only for ^ > 0.8.^max- 

Fig. 3 Surfaces of section (6 = 2) for A = 0.63Acrt (top) and A = 0.76Acrt (bottom). 

Fig.4 The major orbits with winding numbers 5 : 2 (A/Acrt = 0.63) and 3 : 1 (A/Acrt = 
0.76). 

Fig.5 Surfaces of section for 6 = 1 (top) and b — 1/2 (bottom) for A = O.SAcrt- 

Fig.6 Spectrum (b = 2) for A = (a) and for all A (b) as a function of A. (c) and (d) 
are the A = spectra for 6 = 1 and b — 1/2, respectively. 

Fig.7 a.) Contours of SE{N,\) with 100 < < 700 as absissa and < A/Acrt < 1 
as ordinate. Dark areas represent minima. Volume conservation is taken into account in 
Figs. (7) and (8). In b),c),d) and e) the same quantity is displayed for b = 1.5, 6=1,6 = 0.5 
and 6 = 0.58, respectively. 

Fig.8 Shell structure for fixed values of A and 6 = 2. The second derivative of SE{N, A) is 
illustrated versus 

ArV3 fQj, ^ 0, 0.7, 0.63 and 0.76 in (a), (b), (c), and (d), respectively. 

In (e) and (f) the unperturbed (A = 0) structure is displayed for 6 = 2.5 and 6 = 3, 
respectively. 

Fig. 9 Fourier transform of the modulus square of the level density for A = 0, 6 = 2 for 
^/Krt — 0.63 (top) and A/Acrt = 0.76 (bottom). 

Fig.lO Spectrum (6 = 2) for A = extended beyond Acrt = 0.125. Note that the 

distribution of the avoided level crossings occurs largely for A > Acrt- 

Fig. 11 Surfaces of section for the modified Nilsson model. Two situations of accessible 
phase space (shaded) are displayed; left: |A| < uj'^/{AE), right: |A| > uj'^/{AE) 
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